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Abstract. We discuss and review several thermodynamic criteria that have been 
introduced to characterize the thermal stability of a self-correcting quantum memory. 
We first examine the use of symmetry-breaking fields in analyzing the properties of 
self-correcting quantum memories in the thermodynamic limit: we show that the 
thermal expectation values of all logical operators vanish for any stabilizer and any 
subsystem code in any spatial dimension. On the positive side, we generalize the 
results in [R. Alicki et ai, arXiv:0811.0033 ] to obtain a general upper bound on 
the relaxation rate of a quantum memory at nonzero temperature, assuming that the 
quantum memory interacts via a Markovian master equation with a thermal bath. 
This upper bound is applicable to quantum memories based on either stabilizer or 
subsystem codes. 



Submitted to: New J. Phys. Special Issue on "Quantum Information and Many-Body Theory" 



Thermodynamic stability criteria for a quantum memory 



2 



1. Introduction 

Thermal fluctuations pose a serious problem for reliable, passive, information storage 
since any open system eventually reaches a thermal equilibrium state in which all 
encoded information is lost. Fortunately, it has been shown that quantum information 
can be reliably stored for arbitrary long times in, say, a 2D quantum memory [TJ by 
means of active error correction and entropy removal. However, the implementation of 
active error correction implying extensive and fast classical input /output to the quantum 
memory poses a serious (but hopefully not insurmountable) experimental challenge. 

The central idea behind self-correcting classical or quantum memories is to do 
without active error- correct ion and prevent thermalization and build-up of entropy by 
the presence of macroscopic "energy barriers" separating encoded states. 

The idea of such self-correcting quantum memory was first introduced in [1] and 
can be viewed as an extension of the ideas of topological protection developed by Kitaev 
[2]. In [1] it was argued that the 2D surface or toric code (2D Kitaev model) would not 
be a self-correcting memory, but a 4D surface code (4D Kitaev model) generalization 
was presented which would be thermally stable. Unfortunately, three spatial dimensions 
is all the room that the natural world seems to provide. 

It is thus of interest to (1) either come up with models for self-correcting quantum 
memories in 3 or fewer dimensions, or (2) show that low- dimensional quantum physics 
does not allow for such passive stability. The latter possibility would provide evidence 
that genuine quantum phases of nature, such as topological phases, would be confined to 
the domain of finite systems and low temperatures: in the thermodynamic limit thermal 
fluctuations would destroy the quantum order at any nonzero temperature. Such a no- 
go possibility would also lend support to the thought that macroscopic quantum states 
suffer from intrinsic decoherence (see [3] for thoughts in this direction). In this sense 
we believe that the question of thermal stability of a passive quantum memory is one 
of fundamental interest. 

In fact, the thermal stability question can also be viewed as a question about the 
nature of the excitations of the quantum memory model. For 2D topological models 
these excitations are point-like pairs of anyons. If we paraphrase the macroscopic energy 
barrier requirement of [4] in terms of the nature of excitations, it relates to a condition 
that the elementary excitations are extended objects; they are the boundary of a two 
or higher-dimensional surface. 

In [5] the subject of self-correcting quantum memories was brought to the fore. 
Bacon introduced two models, now called the 2D Bacon-Shor code or quantum compass 
model, and the 3D Bacon-Shor code, which are examples of quantum subsystem codes. 
The 3D Bacon-Shor model may or may not be an example of a self-correcting quantum 
memory; it is an open question how to analyze its thermal stability. 

The analysis of the thermal stability of a quantum stabilizer or quantum subsystem 
code model in a thermodynamic sense is the subject of this paper. Let us discuss some 
of the literature on this subject. 
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Necessary criteria for thermal stability of a quantum memory were formulated in 
[I] (see also [6]) in terms of a macroscopic distance of the underlying quantum code (i.e. 
zero temperature topological order) and the presence of macroscopic energy barriers. 
It was shown in that paper that all ID and 2D local stabilizer codes fail to meet 
these criteria. The advantage of this approach is that it allows for very general no- 
go results. A disadvantage is that it does not make contact with any operational or 
thermodynamic expression of thermal stability. In particular, to prove positive results 
on particular quantum memory models, it is necessary to more thoroughly analyze the 
thermodynamics of an open quantum memory. The intuition that underlies the idea 
of a self-correcting quantum memory is that errors of increasing weight should map 
encoded states onto excited states with increasingly higher energy. If the quantum code 
has a macroscopic distance which scales with system-size, then high-weight errors will 
have to happen in order to map one encoded state onto another. But such high-weight 
errors will, if the memory is self-correcting, correspond to high-energy states, hence 
there would be (macroscopic) energy barriers between different encoded states. In the 
second part of our paper, Section [7] we will indeed see that the energy associated with 
high-weight errors corresponding to so-called bad syndromes, will play a crucial role in 
bounding the quantum memory relaxation rate. 

Specific results ruling out the existence of finite temperature topological order for 
e.g. 2D toric code, were obtained in [TJIBIE], using in [8] an interesting finite-temperature 
extension of the topological entanglement entropy. Remarkably, these limitations can be 
overcome by including repulsive long-range interactions with bounded strength. Such 
extensions of the 2D toric code were proposed in [10] and are characterized by a diverging 
relaxation time in the thermodynamic limit. Since the requirement of a macroscopic 
energy barrier between logical states [U |6] is violated in these models, the increase of 
the lifetime with the system size is only polynomial. However, the scaling power is very 
sensitive to the physical features of the thermal bath and becomes especially favorable 
for super-ohmic reservoirs. Such properties needed to be established in [10J by the 
explicit analysis of the non-equilibrium time evolution, instead of being addressed via a 
suitable equilibrium quantity as in the present work. 

In [9] a thermodynamic criterion was presented for the existence of topological order 
at finite temperature. There, it was discussed whether the thermal expectation value 
of logical qubit operators could serve as a stability criterion for a quantum memory 
against thermal fluctuations. Specifically, following the reasoning used in the discussion 
of spontaneous symmetry breaking, a small perturbation (external field) is applied to 
the system which breaks explicitly the symmetry of the Hamiltonian and the state of 
the system. Then, the thermodynamic limit is taken before the external field is taken to 
zero. If the expectation values of the logical operators vanish in this order of limits, then, 
according to the argument given in Ref. [9], the information in the quantum memory 
will be lost after a finite, size-independent relaxation time at any finite temperature. 
This concept was demonstrated explicitly for the Kitaev model in 2D and for some 
generalizations of it to higher dimensions [9]. 
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In this paper we will discuss and analyze this criterion. By making use of elementary 
arguments we show that the same line of reasoning as in Ref. [9] allows one to go well 
beyond these results: in particular, zero thermal averages for the logical operators are 
obtained not only independently of any microscopic details of the code, be it a stabilizer 
code or a subsystem code, but also in any spatial dimension. We will discuss the root 
cause of these problems and discuss possible ways to extend the traditional analysis of 
spontaneous symmetry breaking to detecting a finite temperature quantum order. 

The analysis of thermal stability of a quantum memory within the formalism of 
the thermodynamics of open quantum systems was seriously undertaken in a series of 
papers by Alicki and co-workers [TU [12, [13]. In [12] it was demonstrated that for the 
2D surface code model weakly coupled to a Markovian environment, the relaxation rate 
of any logical state is bounded from below by a constant independent of system size 
[j~2] . This result implies that increasing the system size does not increase the lifetime 
(stability) of the memory, but that the relaxation rate is an intrinsic feature of the 
model. In [13] the authors considered the 4D Kitaev model and explicitly proved that 
the relaxation times were exponentially increasing with system size, hence confirming 
the anticipated thermal stability in the thermodynamic limit. 

In the second part of our paper (Sections [6]TT|) we will present a formal analysis of 
the thermal stability of a quantum memory based on subsystem (stabilizer) codes [14j. 
The difficulty for Hamiltonian models based on subsystem codes (see discussions in [4]) 
is that the Hamiltonian is a sum of non-commuting terms, hence spectral information 
for such systems is not readily available analytically. Using some of the ideas developed 
in [91 [HI [12l [13] we will construct a simple observable whose expectation value on the 
thermal Gibbs state provides an upper bound on the relaxation rate thus determining 
how long quantum information can be stored in a given system. Our formalism will 
be general enough to cover both stabilizer as well as subsystem codes. In addition, we 
can use this formalism to provide a simple bound on the memory relaxation time of 
stabilizer of subsystem code which is not self-correcting, but is 'protected by a gap'. In 
Section [8] we prove that the memory relaxation time scales as n~ l exp(/3A) where n is 
the system size, j3 is the inverse temperature, and A is the spectral gap of the memory 
Hamiltonian. For sufficiently small temperature, e. g. logarithmically scaling with the 
system size, such models may still be of practical interest. 

At the end of the paper we will show that the bound on the relaxation rate only 
depends on an induced temperature-dependent distribution associated with the Abelian 
stabilizer group of the subsystem gauge group. 

2. Stabilizer and subsystem codes 

We assume that the system chosen as the storage medium is represented by an n- 
qubit Hilbert space TC. Let V n = (il, X\, Zi, . . . , X n , Z n ) be the Pauli group on n qubits 
generated by single-qubit Pauli operators and the phase factors ±1, ±1 We envision that 
the quantum data is stored in the degenerate ground states of a quantum Hamiltonian 



Thermodynamic stability criteria for a quantum memory 



5 



acting on the n physical qubits. The Hamiltonian will be associated with a quantum 
stabilizer or subsystem code. 

A subsystem code is determined by its gauge group Q which can be an arbitrary 
subgroup of Q CP n . The set of Pauli operators P <EV n that commute with all elements 
of Q is called the centralizer of Q and is denoted as C{Q). The Abelian group S = QnC(Q) 
is called the stabilizer group of Q. If Q is Abelian, then obviously Q = S up to phase 
factors and we call S a stabilizer code. To preclude S from containing non-trivial phase 
factors one usually adds a requirement —I ^ S in the case of stabilizer codes. If Q is 
non- Abelian, we refer to Q subsystem code. 

Logical operators of a stabilizer code S are elements of C(S) which are not in S. One 
can always choose a set of logical Pauli operators X\, Z\, . . . , X k , Z k G C(S)\S obeying 

2 2 

the usual Pauli commutation relations: X i = Z i = 1 and X; t Zj = (—l) Sl - J ZjXi. Note 
that C(S) = (S,X h Z h ...,X k ,Z k ). 

The code space of a stabilizer code is defined as the common 2 k dimensional +1 
eigenspace of S. It can also be viewed as the ground space of a Hamiltonian acting on 
n qubits: 

m 

H = J2nSi. (i) 

i=l 

Here the are some real negative coupling constants and the operators S% form an 
(over) complete set of generators of S. Note that the definition of a logical operator is 
not unique, since we can multiply any logical operator by an element in S which acts 
trivially on any state in the code space/ground space. Note that the logical operators 
are symmetry operations of H since they commute with all elements Si G S thus each 
energy level of H has a degeneracy 2 k . Therefore, the choice of the ground space as 
coding space, instead of any of the higher energy levels, is somehow arbitrary and other 
forms of encoding might be more useful. An interesting example is the thermal state 
encoding which will be described in Section 16.11 

Bare logical operators of a subsystem code Q are elements of the centralizer C(Q) 
which are not in Q. One can always choose a set of bare logical Pauli operators 
Xi, Zi, . . . , Xk, Z k G C(Q)\Q obeying the usual Pauli commutation relations. Note 
that C{Q) = (S,Xx, Zi, . . . ,X k , Z k ). We can multiply such bare logical operators by 
elements in Q to get so-called dressed logical operators, which act on the gauge qubits, 
in addition to the logical qubits. With the group Q and its local generators Gi we can 

associate a Hamiltonian 

in 

H = ^G t , (2) 
i=i 

where are some real coefficients. Since any Gi commutes with all the bare logical 
operators (Xi, Zi), it follows that H commutes with (Xi, Zi). In addition, H commutes 
with all elements in the Abelian stabilizer group S = Q R C(Q) of Q, hence H is block- 
diagonal in sectors labeled by the quantum numbers (syndromes) of this stabilizer group 
S. Typically, ground states of H are confined to a single syndrome sector. 
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For simplicity we will assume in the remainder of this paper that a single qubit is 
encoded in the quantum memory, i.e. k — 1. 

3. Thermal fragility? 

To get started, let us consider the thermal fragility criterion introduced in [9] and 
apply this to general stabilizer code Hamiltonians, Eq. ([1]). As in Ref. [9], we introduce 
Hh = H — h-S where the additional perturbation is a symmetry-breaking field, designed 
to produce a finite expectation value of the logical operators for the encoded qubit. Here 

s = (x,Y,z). 

For simplicity, we will consider a perturbation h along the ^-direction (this can 
always be assumed by a suitable choice of the logical operators), i. e. h = hn z and 

h • S = tiz. 

We can write the degenerate eigenvectors of if with energy e s as \ s, a) where a = ±1 
is the eigenvalue of Z (the H and the Z operator can be diagonalized simultaneously). 

Clearly, H only acts on the s quantum number (the error syndrome, see Section Hj) of 
the eigenf unctions \s, a), while the logical operators, and in particular the perturbation 
hZ, only acts on the a quantum numbers. As a consequence, the canonical partition 
function Zh = Tr(e _/3jHh ) at temperature 1//3 factorizes 

Z h = ^e-^{s,a\e^\s,a) = (Tre~^) cosh(^) , (3) 

s,a 

where we used '^2 a (s,a\e l3hZ \s,a) = 2 cosh(/3/i), independent of s. 

As was shown in [9], we immediately see that the average value {Z)h is independent 
of the unperturbed Hamiltonian H, and only reflects the finite degeneracy of the energy 
levels 

(Z) h = ^ V- 1 =T — = tanh(/3/i) . (4) 

J2 s ,a e _/3es (s, a\eP hZ \s, a) 

The expectation value in Eq. (jlj) evidently goes to zero at small h. It is clear that 
Eq. (jlj), being independent of the form of H, holds also if the unperturbed Hamiltonian 
refers to a macroscopic system. Therefore, this procedure yields vanishing averages also 
after the thermodynamic limit is taken. If the Hamiltonian H involves n physical qubits, 
we get 

lim lim (Z) h = 0. (5) 

Given that this argument is independent of dimensionality, and thus also holds for the 
4D Kitaev model which is believed to be thermally stable, the result suggests that the 
symmetry-breaking field used in is not strong enough to bias the thermal state 
exp(— f3H h ) towards having a non-zero logical operator expectation value. 

Proof based on the Bogoliubov inequality. We consider next an alternative approach 
based on the Bogoliubov inequality pj3, [T6j and show that we reach the same conclusion 
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as before. This method can then be applied to subsystem codes (see Section [XT]) . We 
start from the well-known Bogoliubov inequality [T5J [16] 

^({A,Ai})([[C,H],tf])>\([C,A})\ 2 , (6) 

where A, C are two arbitrary operators and H is the system Hamiltonian, with the 
assumption that all expectation values exist. Here we use the convention {A, B} = 
AB + BA. We then set A = X and C = Y. Clearly, {X,X} = 2 and the right-hand- 
side of the Bogoliubov inequality (jSJ) gives 4\(Z) h \ 2 . Therefore we get 

P([[Y,H-hZ},Y]) h >A\(Z) h \ 2 , (7) 

where the Y operator commutes with H. Using the commutation relations for the logical 
operators, we obtain A(3h{Z) h > 4\(Z) h \ 2 > 0. For (Z) h strictly positive (otherwise we 
are done) we can divide by {Z)h, and then, by taking the thermodynamic limit on both 
sides of the resulting inequality, we eventually get 

lim (Z) h < (3h . (8) 

n— >oo 

At any finite temperature, we thus obtain that the thermal expectation value of the 
logical operator vanishes when h — *■ 0. 



3.1. Subsystem codes 

Let us use the Bogoliubov inequality to argue about the thermal fragility criterion 
for subsystem codes. The bare logical operators of the encoded qubit S = (X, Y, Z) 
commute with all Gi G Q, hence with H in Eq. (j2J). We can choose the symmetry- 
breaking Hamiltonian as 

H h = H — hZG , (9) 

for some choice of G which dresses the bare logical operator Z. Let us thus consider the 
thermal expectation value of (ZG')h where G' does not need to be the same as G. 

We use the Bogoliubov inequality, with A = XG' and C = Y. Since G' commutes 
with X and G'G'^ = 1 (valid for every Pauli operator), one obtains {XG', G'^X} = 2. 
This gives 

f3([[Y,H - hZG],Y]) h > 4\(ZG') h \ 2 , (10) 

and since Y commutes with H and G, we can easily compute the left-hand side, to 
obtain 

Ph(ZG) h > \(ZG') h \ 2 . (11) 

We now notice that ZG is a Pauli operator and thus has eigenvalues ±1. Hence, the 
thermal expectation value on the left side is always less then 1 independently of any 
details of H, which gives 

lim \(ZG') h \ 2 < (3h , (12) 

n— >oo 

for any choice of G and G'. Therefore, the same considerations valid for the stabilizer 
codes can be repeated in this case and we again conclude that the thermal expectation 
value of any logical operator vanishes at any finite temperature for vanishing field h. 
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4. Error correction 

Let us pause for a moment and discuss our somewhat naive-looking approach. It seems 
that there are at least two issues at stake here. Let us assume that by choosing the right 
symmetry-breaking field, we are able to concentrate the weight of exp(—/3Hh) around a 
logical |0). Consider this +1 eigenstate |0) of the Z logical operator and a state with a 
single qubit error, E\0), such that E anti-commutes with Z. Obviously, if at equilibrium 
the system is in a statistical mixture of |0) and E\0) with equal probability, one has 
(Z) = 0. However, the information in the memory is still preserved as long as we 
correct for errors such as E when we determine what logical state has been stored. For 
a generic stabilizer code, the probability of the E\Q) states is small at low temperature 
(below the gap), but the statistical weight of all correctable errors might be very large in 
the thermodynamic limit. Therefore, (Z) does not represent a meaningful stable order 
parameter for this problem; the value of Z has to be modified depending on the error 
syndrome. For this reason, the authors of Ref. [13J consider so-called error- corrected 
logical operators^. Let us properly define these for stabilizer and subsystem codes. 

For stabilizer codes, error correction consists of measuring the ±1 eigenvalues of 
the stabilizer generators; these sets of eigenvalues form the error syndrome. The error 
syndrome is used as input to a classical decoding algorithm which determines which 
errors have most likely taken place. For subsystem codes, error correction may proceed 
by measuring the eigenvalues of the local generators Gj. Since the operators Gi do 
not commute, these eigenvalues cannot be simultaneously measured, nonetheless these 
(random) values of the generators of Q will fix the eigenvalues of the stabilizer group S. 
These eigenvalues of the stabilizer group S again form the error syndrome. 

More precisely, any error E £ V n determines a syndrome se '■ <S — > Z2 such that 



We can assume that there is some deterministic decoding algorithm which assigns a 
correcting Pauli operator C(s) £ V n to every syndrome s. An error E £ V n is correctable 
iff C(se) coincides with £ up to a gauge operator, that is, EC(se) £ G- 

We can define a subspace projector P s associated with every syndrome (quantum 
number) s. Let Pq be the projector onto the 5-invariant code space in which Si\ip) = 
for all i — 1, . . . ,p. (By abuse of notations let us assume from now on that —I ^ S.) For 
any syndrome s we can define P s = EPqE^ where E £ V n is any error with syndrome s 
(note that the projector P s does not depend on the choice of such E). Clearly J2 S — I- 
We define an error- correcting transformation for observables on TC as 



Note that $ ec (J) = I, so the adjoint transformation $* c acting on states is a trace- 
preserving completely-positive (TPCP) map. Following [13] we can define the error- 

| In [T3] these are called dressed logical operators, but we prefer to reserve the notion of 'dressing' for 
the multiplication of bare logical operators with elements of the gauge group Q. 



EQ = (-l)^W) QE for all Q £ S. 




(13) 



s 
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corrected logical operators as 

Z ec = $ ec (Z), X ec = $ ec (X). (14) 

for a pair of bare anti-commuting logical operators (X, Z). Note that Z ec ,X ec are not 
necessarily Pauli operators. However, it is not hard to show that the error-corrected 
logical operators obey the relations Z ec = X ec = I and Z ec X ec = —X ec Z ec . We can 
understand this by defining coefficients X z (s), X x (s) G {+1, —1} such that 

C{sYXC{s) = X x (s)X, C{sy~ZC{s) = X z {s)Z. (15) 

Any syndrome projector P s belongs to the algebra generated by S and thus commutes 
with Z, X. It follows that 

Z ec = ZD Z , where D Z = J2 Ps, (16) 



X ec = XD X , where D x = 2J X x (s) P s . 



The commutation relations for Z ec ,X ec follow directly from Eq. ( |T6i) . Note also that 
the error-corrected logical operators commute with all elements in Q. 

We can immediately check whether the use of error-corrected logical operators 
would change the analysis of the thermal expectation values. As observable, we choose, 
say, Z ec G' for some G' whereas for the symmetry-breaking field we choose some Z ec G. 
Using the properties of Z ec stated above, we can repeat the proof of the previous 
subsection to obtain again a vanishing expectation value 

lim lim (Z ec ) h = 0. (17) 

h— »0 n^oo 

5. Analogy with the 2D Ising model: choice of symmetry-breaking field 

We emphasize that the conclusions above are valid for arbitrary dimensions of any 
stabilizer or subsystem code. Although (or since) the argument is so universal it also 
appears to be exceedingly oversimplified. In the previous section, we have discussed 
the necessity to choose a stable logical observable which includes the process of error 
correction. Let us now more closely examine the choice for the symmetry-breaking field. 

Although the thermal fragility criterion is patterned along the lines of standard 
symmetry-breaking arguments, it is only so on a formal level. It is instructive to compare 
the argument of Ref. [9] with the standard example of spontaneous symmetry-breaking 
in the 2D Ising model [19] (see e. g. [20]): 

where i, j label the 2D sites of the full lattice A, the first sum is over pairs of 
nearest neighbor sites, and an external magnetic field b is included. For the 2D Ising 
model one obtains at low temperature 

lim lim (Zj) b ^ , (19) 

b— >0 n— >oo 
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at every lattice site j, where the expectation value above is taken with respect to the 
Hamiltonian (IT81) . This appearance of a symmetry-breaking order should be contrasted 
with the lack of such order in the ID Ising model which has T c = 0. 

Notice that, although the 2D Ising model does not display topological order, it does 
define a proper stabilizer code with logical operators X = Ylie\Xi anc ^ % = Zj, where 
j is a fixed (arbitrary) site in the lattice^. 

The arguments discussed in the previous sections consider a perturbation hZ = hZj 
which leads to (Z)h = (Zj)h — > and does not show that the value of the ^-polarization 
is robust. In fact, the field h only acts on a single site, whereas in the standard case the 
symmetry breaking field b acts on all sites of the lattice A simultaneously, see Eq. ( 1181) . 
The reason for the failure of the stability criterion appears thus to be that the chosen 
symmetry-breaking perturbation is not extensive. Although for topological memories 
the support of a logical operator Z (i.e., the number of physical spins on which the 
operator acts nontrivially) becomes larger with the size of the system, the perturbation 
hZ is bounded in norm by h and becomes irrelevant in the thermodynamic limit. 

The analogy with the 2D Ising model suggests that the symmetry-breaking field 
should be chosen as a sum over different incarnations of a logical operator, i.e. we 
can multiply a logical Z by elements of the stabilizer code S and obtain an extensive 
operator. It may be possible to salvage this symmetry-breaking route to getting 
a quantum order parameter, but of course any construction should ultimately be 
motivated operationally. This is the reason that we now switch to explicitly deriving a 
memory relaxation rate. 



6. Relaxation rate for general quantum memory Hamiltonians 

The goal of this section is provide a criterion for thermal stability for a large class of 
quantum systems that can be described by subsystem codes [H] . This is a generalization 
of the work in Ref. [13] in which the thermal stability of the 4D Kitaev model was 
analyzed by considering the dynamics of the quantum memory in contact with a thermal 
bath. 

Let TC be the Hilbert space describing the system chosen as a storage media and A 
be the algebra of operators acting on TC. The following definition will play an important 
role in this section. 

Definition 1 Let O G A be an observable and let P be a projector onto some subspace 
ofTC which is invariant under O, that is, PO = OP. We shall say that the observable 
O is protected from a set of errors £ C A on a subspace P iff 

[E, O]P = for all EeS. (20) 

§ Of course, the other two expectation values (X) and (Y) are vanishing in the appropriate 
thermodynamic limit. This stabilizer code does not provide a good quantum memory since the distance 
of the code is 1 independent of lattice size. 
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(Here and below we use the notation P both for a subspace and the corresponding 
projector.) Consider as example the case when S includes all single-qubit Pauli 
operators. Suppose O \ip) = A|^) for some \if)) G P. Then Eq. ( 1201) implies that 
OE = XE for all E G S, that is, a single-qubit error cannot change the eigenvalue 
of O for any eigenvector that belongs to P. Quantum error correcting codes provide 
a systematic way of constructing observables protected from low-weight errors on a 
code-subspace, see below. 

Suppose for simplicity that our goal is to encode a single qubit. We shall need a 
pair of observables X, Z G A obeying the canonical commutation rules of the Pauli 
operators, 

X 2 = I, Z 2 = I, XZ = -ZX. (21) 

In the following we shall refer to X and Z obeying Eq. (l2Tj) as Pauli-like observables. 
(Note that Pauli-like observables need not to be single-qubit Pauli operators or tensor 
products of Pauli operators.) 

Assume that the system evolves according to a Markovian master equation 

p=-i[H,p]+£{p), (22) 

where C : ^4. — ^ ^4 is the Lindblad operator defined by 

£ = ^£ £ a (p) = S aP Si-±{p,SiS a }. (23) 

a 

The operators S a will be referred to as quantum jump operators. For any Lindblad 
operator C, let Sc C A be the set of all quantum jump operators involved in C. 
Integrating Eq. (1221) one arrives at 

pit) = $*(p(0)), $ t = exp (-it[H, ■] + tC). (24) 

We shall measure the strength of £ using the norm 

||£||i = max ||£(F)||i subject to ||F||i < 1. (25) 

Here the maximization is over all self-adjoint operators F = acting on the system 
Hilbert space and ||F||i is the trace norm of F, i. e. = TrvFFt. Note that \\F\\i 

is distinct from the spectral norm \\F\\. 

The following theorem is the main result of this section. 

Theorem 1 Let C be an arbitrary Lindblad operator with a set of quantum jump 
operators £c such that the Gibbs state pp ~ exp {—(3H) is the fixed point of C, C{pp) = 0. 
Suppose one can choose Pauli-like observables X, Z G A that are protected from the set of 
errors Sc on some subspace P. Suppose also that X, Z , and P commute with the system 
Hamiltonian H . Then there exist TPCP encoding and decoding maps $j n : L(C 2 ) — > A 
and $ OU ( : A — > L(C 2 ) such that 

||$t o $ in (7j) - $ in (ri)\\i < St \\C\U Tr(/ - P) Pf3 (26) 

and 

®out o <$> in (r}) = 7] (27) 
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for all one-qubit states r] and for allt > 0. 

Note that the right-hand side of Eq. (126]) provides an upper bound on the precision up to 
which the decoded state Q ou toQ t oQ in (7]) approximates the initial state rj. Thus assuming 
that the system consists of n qubits and that the norm of the Lindblad operator grows 
at most as poly{n) we can store a qubit reliably for a time of order 

c qmem ) 

(28) 

where 

eqmcm = Tr(J - P) Pf3 . (29) 

We shall refer to r qmem as the storage time and to the quantity e qmem as the relaxation 
rate. One can envision two scenarios when the bound Eq. (128]) on the storage time can 
be useful: (i) the relaxation rate e qmem is exponentially small as a function of n, that is, 
eqmem < ex P ( — nl ) for some 7 > 0; (ii) the relaxation rate e qmem is only polynomially 
small but the degree is sufficiently large, such that r qme m grows fast with n. The first 
scenario can be realized for systems featuring a macroscopic (growing as n 1 ) energy 
barrier surrounding the states orthogonal to the protected subspace P. The 4D toric 
code model analyzed in [13] provides an example of such a system. The second scenario 
could be realized if the energy barrier grows only logarithmically as a function of n 
as in [17] HE]. In this case the exponent is controlled by the temperature, that is, 
eqmcm < n^" 1 ^ = e~ ^ log ^ for some 7 > 0. If the temperature is smaller than a critical 
value, the relaxation rate e qmem decays sufficiently fast to yield a storage time r qmem 
increasing with n. A polynomial increase of the storage time is also obtained in [10] at 
any temperature, from the logarithmic divergence of a self-consistent gap. It is tempting 
to conjecture that such system may exist in lower spatial dimensions. 

The proof of Theorem [1] involves two ingredients: (i) constructing the encoding 
and decoding maps (see Section 16.11) . and (ii) proving that the encoded states are 
approximate fixed points of the Lindblad operator (see Section I6T21) . Our construction 
of encoding and decoding maps is identical to the one used by Alicki et al. in [Hi H3] . 
It is described in Section |6J] which can be regarded as an overview of Section IA in [13J. 
The second part of the proof is presented in Section 16.21 Our approach here is quite 
different from the one taken in [13]. It yields a much simpler proof and requires less 
assumptions about the Lindblad operator compared to [13] (for instance, we don't need 
the detailed balance condition). 

Following pJ!}[12|,[T3] we ca n specialize Theorem[T]to the Markovian master equation 
due to Davies [21] which describes the dynamics induced by a weak coupling between 
the system and a thermal bath. It involves a coupling Hamiltonian 

K 

H int = ^A k ®B k , (30) 
fe=i 

where A k are some local few-qubit operators acting on the system and the operators B k 
act on the bath. 
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It was shown by Davies [2TJ that in the weak-coupling limit the system evolves 
according to the Markovian master equation Eq. (1221) where the Lindblad operator is 
defined as 

C ^ = E E ( a ^p a L - 1 4,oA, «}) • (si) 

k LU ^ ' 

Here are the Fourier components of Ak(t) = e lHt A^e~ %m , that is, 
A fc (t) = ^A fciW e- iwt . 

id 

One can think about u as the part of Ak transferring energy to from the system to the 
bath. The bath temperature enters into the equation only through the function h(k, to) 
which has to obey the detailed balance condition, 

h{k,-u) = e~ /3ui h{k,u}). (32) 

The coefficient h(k, to) is defined as the Fourier transform of the autocorrelation function 
of Bk with respect to the bath state. One can regard h(k, to) as a probability (per unit 
of time) of quantum jumps induced by the coupling operator Ak which transfer energy 
to from the system to the thermal bath. The detailed balance condition guarantees that 
the Gibbs state pp is a fixed point of C 

It is important to discuss how the quantum jump operators depend on the 
original coupling operators Ak- 

For stabilizer code Hamiltonians as in Eq. (pQ) the time- dependent operator Ak{t) = 
exp(iHt)Akexp(—iHt) acts only on a few qubits since all the terms in H pairwise 
commute and thus Ak(t) = exp(iH't)Ak exp(-iH't) where H' includes only those terms 
of H that act on the same qubits as Ak- Note that H' has only a few Bohr frequencies 
since it acts only on a few qubits. It means that any quantum jump operator Ak^ in the 
Davies master equation acts only on a few qubits and the total number of the quantum 
jump operators is roughly the same as the number of the coupling operators Ak- 

This issue is more subtle for subsystem codes, since Ak(t) may be a highly non- 
local operator for long times t and the number of Bohr frequencies may be exponentially 
large. However, it is also clear that the non- locality of Ak(t) is only due to multiplying 
it with non-local elements in the gauge group Q. Hence Ak(t) remains local modulo 
gauge group transformations. 

Let us specialize the Theorem [1] to the Davies master equation, see Eqs. (I22ll3ip . 
The condition that the observables X and Z are protected from all quantum jump 
operators in Sc might seem too demanding since the operators A^ u may be highly non- 
local, see the remark above. Fortunately, it is sufficient to require that X and Z are 
protected from a set of errors S int = {A k } including all coupling operators A k . Indeed, 
since, by assumption, H commutes with P and X, Z, the condition [X, Ak] P = 
implies [X, Ak iLU ] P = for any frequency to. (The same remark applies to Z.) 

Next we need an upper bound on the norm of the Davies generator C, see Eq. ( 13~T|) . 
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Proposition 1 Assuming that \\A}~\\ < 1 for all k,u one has 

\\C\U < 2Kh max , (33) 

where h max = maxj;^ \h(k,uj)\ and K is the total number of terms in the interaction 
Hamiltonian Eq. [Sty) . 

Proof. Indeed, let F = be an operator such that \\F\\i < 1 and ||£||i = ||£(F)||i, 
see Eq. f l25l) . Fix some k and let A = Ak, A u = Ak }0J , and h{uj) = h{k, uj). Let us bound 
the trace norm of a single term 

c k (F) = Yl K") A « FA l - {A « A «> F} - (34) 



Note that 



and 



4A = T lim -L dt A(t)^A(t) (35) 



J^A^FAl =lim ± £ dtA(t)FA(ty. (36) 

uj ~T 

Using the bound ||AB||i < ||A|| ■ valid for any operators A,B we get 

11^(1/2)^(0;) F}\U < 11^^)4^11 ■ \\F\U (37) 

UJ UJ 

^ hmax || ^ A^A^ || • || .f 1 ||i ^ h ma x 1 1 F 1 1 1 ■ 
uj 

Here the second line used Eq. (135"!) . convexity of the norm, and the fact that ||A(£)|| = 
Pll < !• 

Let F = F + — F_ be the decomposition of F into positive and negative parts, that 
is, F± > and ||F||i = Tr F_ + Tr F + = ||F_||i + Then 

H^HA^H^ || ^/ i (o;)^F_4|| 1 + 11^^)^+411! (38) 

a; a; a; 

5: ^mosll ^ ] A^F_y4jj| 1 + h max || ^ ^ A^i^A^H i 

a; uj 

^ hmax 1 1 F_ 1 1 1 "I - h max 1 1 1 1 x h max 1 1 F 1 1 x . 

Here the last line used Eq. ( 136]) . convexity of the norm, and inequality ||A(£)i<±A(t)t||i < 
II -Ft ||i- Combining Eqs. (I37II38I) we arrive to ||£fc(F)||i < 2/i moa; ||F||i which leads to 
Eq. (j33J. 

□ 

To conclude, Theorem[T]can be specialized to the Davies master equation as follows. 
Suppose the system interacts with a thermal bath at the inverse temperature (3 via a 
Hamiltonian H int = Ylk=i Ak ® wnere ll^fell — 1 an d B k are normalized via the 
condition h(k,u) < h max . Suppose one can choose Pauli-like observables X,Z that are 
protected from any coupling operator A k on some subspace P. Suppose that X, Z, and 
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P commute with the system Hamiltonian H. Then Theorem [T] implies that a qubit can 
be stored in the system reliably for a time r qmem ~ (M ma3 .e qmem ) _1 , where 



-qmcm = Tr(I - P)p p . (39) 

Note that K will be 0{n) for local couplings A^. We will discuss how to evaluate the 
relaxation rate e qmem in more detail in Section [7J 

6.1. Proof of Theorem^ part I 

Let us start from defining the encoding and decoding maps $j n and $out- Let Aq C A 
be the algebra generated by J, X, Z, and K = iXZ. For any algebra A let us define 
the center of ^4 as 

= {A e A : AB = BA for all B e A}. 

Clearly Z{Aq) = C ■ I, that is, Aq has trivial center. 

For any finite-dimensional Hilbert space let h(H) be the algebra of linear operators 
acting on H. We shall use the following fact (see for instance Theorem 5 in [22], or a 
book [23J): 

Fact 1: Let Aq C ~L(H) be any algebra such that (i) Aq contains the identity operator; 
(ii) Aq is closed under hermitian conjugation; (Hi) Aq has a trivial center. Then there 
exists a (virtual) tensor product structure H = TCq ® Ha such that 

Aq = h(H Q ) ® I A . (40) 

It implies that there is a decomposition H = Hq <S> TLa such that TCq describes a 
qubit Q and the operators X, Y, Z are the Pauli operators acting on TIq, that is, 

X = X Q ®I A , Y = Y q ®I a , Z = Z q ®I a . (41) 

By assumption, the system's Hamiltonian H commutes with X, Y, Z. Therefore H acts 
trivially on Hq and thus there exists Ha € L,(Ha) such that 

H = I Q ® H a . (42) 

Note that Trexp (—f3H) = 2Trexp (—(3 Ha). Thus the Gibbs state pp can be written 

as 

2 Tr exp (— p f/^J 

Define the encoding map $j n : L(C 2 ) — > A as 

$m(?7) = V®VA- (44) 

Using Eqs. (14114211431) one gets 

$ m (J) = 2 P/3 , $ m (Q) = 2Qp^ = 2 P/3 g for any Q e {X, Y, Z}. (45) 

Define the decoding map <& out : A — > L(C 2 ) formally as the partial trace over the 
subsystem Ha, 

^out(p) — Tr^ p. (46) 
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Clearly, & out o $ in is the identity map. 

To demonstrate this formalism, let us explain how the encoding map $j n is 
constructed for the special case of stabilizer (subsystem) codes. Imagine that one needs 
to store a single qubit state r\ = ~(I + 7 • S) with S = (X, Y, Z). We encode into the 
thermal state $i„(r]) = 2ppr] ec where r] ec = ~(I + 7 • S ec ) with the error-corrected logical 
operators S ec = (X ec , Y ec , Z ec ). Note that rj ec commutes with pp. 

The central idea underlying the encoding into the thermal state is that &i n (r]) is 
the same as the stationary state pp satisfying Qtipp) = p@ on Ha- Thus we can expect 
that if thermal fluctuations do not build up to logical errors, the state $ t o ^ in (rj) would 
remain close to the initial state $i n (v)- 

Note that $« n is quite different from the standard encoding into the the ground 
state subspace, for which the requirement of a Hamiltonian with finite excitation gap 
appears most natural. Instead, the stability criterion of Theorem [1] using the encoding in 
a thermal state does not explicitly involve the spectral gap. This is an interesting point, 
since it has become clear now that the presence of a gap does not imply robustness of 
topological protection. On the other hand, it might be possible to obtain a self-correcting 
quantum memory for a Hamiltonian with vanishing gap at large n. 

6. 2. Proof of Theorem [IJ- part II 

Let p = ^inijf) be any encoded state. Using Eq. (05]) one check that p can be represented 
as 

p = Opp = p p O, OeAq, \\0\\<2. (47) 
Consider a family of states 

p(t) = $ t (p), $ t = exp(-it[H r ]+£t), t>0. (48) 
Taking into account that H commutes with p we can represent the derivative p as 

p(s) = $ a (£(p)), s > 0. (49) 
Using the fact that ||$(A)||i < ||A||i for any TPCP map $ and any operator A we get 

||/K*)i|i<||£(>)||i, s>0. (50) 

Therefore 

I|p(*)-P(0)lli = ll f dap{s)\\i<t\\^{p)\\i- (51) 

Thus we have to prove an upper bound on the norm of C(p) = C(Opp). Inserting twice 
the decomposition P + P 1 - = I we get 

C{Op p ) = Li + L 2 + L 3 + L 4 , where (52) 
L x = P£(OPpp), L 2 = P ± £(OPpp), 

and 

L 3 = P£(OP ± pp), U = P L £(OP ± pp). 
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Using the identity ||AB||i < ||A|| • \\B\\i valid for any operators A, P, taking into account 
that ||0|| < 2 and using Eq. ff25l) . one easily gets 



\\L 3 \\ 1 ,\\L 4 \\ 1 <2\\£\\ 1 TrP ± p p . (53) 

We shall bound the norm of L x and L 2 using the fact that C(p@) = 0. Indeed, using the 
assumption that [S a , O] P = and P[Sl,0] = for all a one can rewrite L\ as 

In = OPC(Ppfs) = -OPC(P^pp). (54) 

It follows that 

||Li||i < 2||£|| 1 TrP ± p/3- (55) 
Using P[0, Si] = 0, PP 1 - = 0, and C{pp) = we can rewrite L 2 as 

L 2 = P ± £(Ppp)0 = -P ± C(P ± p^)0 (56) 

and thus 

||L 2 ||i <2||£|| 1 TrP ± p /3 . (57) 
Combining Eqs. (I53|55ll57p we arrive at 

WCiOp^W^SWCh TrPV (58) 
Plugging it into Eq. (JED) we get \\p(t) - p\\ x < 8||£||i Tr(J - P)p p . 

7. Relaxation rate for subsystem code Hamiltonians 

In this section we explain how to construct the protected Pauli-like observables and 
the subspace P involved in Theorem [1] using the formalism of subsystem codes. Let 
Q Q V n be the gauge group of some subsystem code encoding one qubit into n qubits. 
Assume that the system's Hamiltonian H is defined as in Eq. (j2j), so that if is a linear 
combination of gauge operators. Suppose we seek protection from some set of elementary 
errors 8. Assume without loss of generality that all elements of £ are Pauli operators, 
that is, £ C V n . For example, £ may include all Pauli operators that appear in the 
decomposition of the operators Ak coupling the system and the bath, see Eq. ( 1301) . In 
this case any elementary error acts only on a few qubits. 

Let us start from defining a notion of goodness of syndromes relative to the set of 
elementary error £. We will say that a syndrome s, see the definitions in Section HI is 
good iff 

C(s + se)EC(s) G g for all Ee£. (59) 

Remember that C(s) G V n is the correcting Pauli operator for a given syndrome s 
which is determined by some deterministic error correction algorithm. To highlight the 
intuition behind the definition of good syndromes, let us assume that the syndrome s 
has been caused by some pre-existing error E'. If the error E' is correctable then we 
have C(s)E' G Q and thus EC(s) coincides with EE' up to a gauge operator in Q. Note 
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that EE' has syndrome s + se- Eq. (15T)j) says that C(s + se)EE' G (/, that is, the error 
is also correctable for all elementary errors E G £. 

We would like to point out that in the theory of quantum fault-tolerance and error 
correction, very similar notions are used to determine the correctness of an encoded 
logical gate (called a rectangle), see e.g. in the discussion at the bottom of page 12 
in Ref. [24J. The correctness of the encoded logical gate depends on its incoming pre- 
existing syndrome in combination with new errors which occur during the execution of 
the encoded gate. In our definition of goodness, no gate happens, but we allow for any 
elementary error E and determine whether the pre-existing syndrome in combination 
with the new error leads to making a good inference about the total error. 

Thus one can keep adding more and more elementary errors as long as the observed 
syndromes are good. On the other hand, if the observed syndrome becomes bad (that 
is, not good), it means that one has already reached the limits of the error correcting 
capabilities of the code and the next elementary error can potentially destroy the 
encoded information. In this case the operator C(s + se)EC(s) G C(S)\Q becomes 
a non-trivial logical operator. 

Now, in order to apply Theorem [U we pick some bare logical Pauli operators 
X, Z G C{G)\Q, see Section [21 and choose the Pauli-like observables X and Z as the 
error-corrected logical operators X ec and Z ec defined in Section HI Below we shall prove 
that X ec and Z ec are protected from the set of elementary errors £ on the subspace P 
spanned by good syndromes, that is, 

p = E p » ( 6 °) 

good s 

Note that by construction X ec , Z ec , and P commute with any Hamiltonian H made up 
from gauge operators, see Section @J while P is an invariant subspace of X ec and Z ec , 
see Section HI so all the conditions of Theorem [1] are met. 

Lemma 1 The observable Z ec is protected from all elementary errors on the subspace 
spanned by good syndromes. 

Proof. Indeed, let s be any good syndrome and E G £ be any elementary error. It 
suffices to prove that 

[E,Z ec }P s = 0. (61) 

Let t = s + se and A = EC(s)C(t). Goodness of s implies that A G Q. Since we have 
chosen Z G C(Q), it implies AZ = ZA. Taking into account Eq. ( fl5i) we get 

ZE = EZX z (s)X z (t). (62) 

Thus we see that 

[E, ~Z ec ] Ps = Pt[E,Z ec ] P s = P t (X z (s)EZ - X z (t)ZE) P s = 

where we have used Eqs. ffTo^ . ffTB]) and (16"2"|) . 

Recall that the syndrome subspaces P s are well-defined only if one fixes the subspace 
P associated with the trivial syndrome. A natural choice of P is dictated by the ground 
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state of the system's Hamiltonian H = J2i r iGi- Typically the degeneracy of ground 
subspace of H is the minimal degeneracy consistent with the symmetry of H. In our 
case the non-abelian symmetries of H include the bare logical operators X and Z, 
so we should expect the ground state to have degeneracy 2. In this case the ground 
state determines a particular set of quantum numbers (irreducible representation) of 
the stabilizer group S = Q H C{Q). We can choose the trivial syndrome subspace Pq as 
the subspace spanned by all states that have the same quantum numbers as the ground 
state. Equivalently, P includes all states that can be obtained from the ground state 
by applying gauge operators and logical operators. 

The question of whether a particular family of subsystem codes is suitable for 
building a good quantum memory can now be reduced to bounding the relaxation rate 
defined in Eq. ([29]): 



e q mcm = 1 ~ TrppP = ^ ds f\ where Z s = Tr P s exp {-/3H), (63) 

as a function of n. Recall that all terms in H commute with S so that all syndrome 
subspaces P s are invariant under H. It is clearly desirable to have a Gibbs state 
exp(— (3H) with support mostly concentrated on the good syndromes. Note that low- 
weight correctable errors which in addition remain correctable if any single additional 
error occurs, will have syndromes which are good. Hence the Gibbs state should be 
concentrated on the ground space and excited states that can be created from the 
ground state by these correctable errors. This type of property is what has been shown 
for the 4D Kitaev model in [T5] . Bounding the relaxation rate for subsystem codes is 
a hard task, since it depends on the full partition function of the model. The following 
simple observations might be helpful for obtaining upper bound on e qmem - 

Consider the partition function Z s associated with some syndrome-sector P s . Let 
E G V n be some error with the syndrome s. Note that P s = EP E. It implies 

Z s = Tr EP E exp (-/3H) = Tr P exp (—(3EHE) 

= Tr P exp (-/3H + (3 He) < Tr P exp (-(3H) exp (J3H E ) , (64) 

where 

H E = 2 r ^ ( 65 ) 

i:G,E=-EGi 

is a sum of all terms in H that anticommute with the error E. The second line in 
Eq. (164jl follows from the Golden-Thompson inequality. It implies 

Cqmem < (Pq exp {(3H E { S )))(1 (66) 



bad s 



where is the average over the thermal Gibbs state and E(s) is some fixed error 
causing syndrome s. 

Let us now express our intuitive understanding of under what circumstances the 
bound in Eq. fl66|) could give rise to self-correction. For a self- correcting quantum 
memory, we expect that bad syndromes correspond to errors which anti-commute 
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with a macroscopic number of terms in H, hence He is a sum over a macroscopic 
number of terms, say l(n). This type of requirement has been expressed in Ref. [I]. If 
this requirement is fulfilled, one can imagine that in 'sufficiently high dimensions', it 
is possible to use a mean-field approximation and approximate (P exp (/3He( s )))/3 by 
exp ((3{PqHe(s))p)- It is of course important that fluctuations around such mean-field 
approximation die off sufficiently fast. Now, if one can upper bound each individual 
term in {PqH e(s)) p by some constant — c, then, because there are a macroscopic number 
of terms in He( s ), the factor (P exp (/3He(s)))/3 would scale like exp (— f30{l{n))). The 
sum of the bad syndromes will multiply this exponential decay by some factor upper 
bounded by 2 rank ( 5 ) where rank(iS) is the minimal generating set of S. If l(n) grows 
at least as fast with n as rank(iS), then e qmem would be exponentially decaying. It is 
clear that many things have to 'go right' in order for self-correction to be feasible, in 
particular it is not clear whether three-dimensions would be sufficient for mean-field 
type approximations with sufficiently small corrections. 

Another simple observation shows that one may be able to make headway in 
computing the partition functions Z s by making use of the underlying symmetry. Recall 
that P s is the projector associated with some irreducible representation of the stabilizer 
group S = Q H C(Q). Therefore we can write P s as a linear combination of elements of 
S, 

P s = J2°s(Q)Q> <Q) £{+!,-!}• (67) 
QeS 

Expanding the exponent exp {—(3H) in powers of f3 k we will only get non-zero 
contributions to Z s from terms in H k that are elements of the stabilizer group S. We 
will leave a detailed analysis of the memory relaxation rate for a particular subsystem 
code to a future paper. 

8. Relaxation Rate For Gapped Hamiltonians 

It is worth emphasizing that the lower bound on the storage time obtained in Sections [6H71 
applies to both gapped and gapless memory Hamiltonians. It is natural to ask whether 
a stronger bound can be obtained if the memory Hamiltonian H has a constant spectral 
gap A separating the ground state from excited states while the bath temperature T is 
small compared to A. In this low-temperature regime the rate of all processes exciting 
the system from the ground state is suppressed by the Boltzmann factor e _/3A and thus 
one should expect that the storage time scales as 

— *— < 0(n) e~P A . (68) 

Tqmem 

Below we prove that this is indeed the case assuming that the interaction with the 
thermal bath can be described by the Davies equation, see Eqs. f)22|3ip . In contrast to 
our main result (see Theorem [1]), the proof of Eq. ( 1681 will rely on the detailed balance 
condition Eq. (1321) . Our analysis will use the encoding and decoding maps § out 
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defined in Section IBTTl where we set (3 = oo. In other words, we encode information into 
the ground state of H rather than the thermal state. Accordingly, any encoded state 
p = <&in(ji) can be represented as 

P = Poo = Poo O, O G Aq, \\0\\<2. (69) 

The observable O must be protected from all coupling operators A k , see Eq. (1301) . on 
the ground-state subspace P , that is, one must have [O, A k ] P = for all k. As was 
mentioned in Section [6], this is equivalent to the condition 

[O, A k> J\ P = for all k, for all u. (70) 

Using Eq. ( I3T1) it suffices to prove that 

||£(p)||i< 0(n)e-^ (71) 

where C is the Davies generator defined in Eq. ( 13TI) . Let Afc ja; be any quantum jump 
operator from Eq. (I3TI) . Since transfers energy uj from the system to the bath, we 
have 

A k , UJ p oo = forcu>0. (72) 

It follows that C(p) contains only the terms with non-positive Bohr frequencies. We 
claim that the zero-frequency term also does not contribute to C{p). Indeed, since Ak,o 
commutes with H, we get 

Ak,oPo = PoAk t o = PoAkfiPo- (73) 
In addition, since A k is hermitian, we have 

4,o = 4fc,o. (74) 
Combining Eqs. (I70|73ll74p one can easily check that 

Ak,opA{ --{p,Al A kt o} = 0. (75) 
Since any negative Bohr frequency is separated from by the gap A, we get 

£ ^) = E E Hk,u) U ki „pAl„-±{p,Al„A k ^}Y (76) 

k uj<-A V - / 

The detailed balance condition Eq. ( |32l) implies that for any uj < —A one has 

\h(k,u)\ <e-P A \h(k,-u)\ <e~^h max . (77) 

Let CJ be a superoperator obtained from C by setting h(k, uj) = for u> > —A. Using 
Proposition [1] and Eq. (1771) one infers that ||£'||i < 2Kh max e~ l3A for some K = 0(n). 
It means that 

||£(p)||i = ||£'(p)||i < ||£'||i < 2Kh max e~ f3A = 0(n)e^ A . (78) 

It proves Eq. (17X1) completing the proof of Eq. fl68l) . 

It is worth pointing out that the condition [O, A k ] Pq = used in the above analysis 
is satisfied whenever the coupling operators {A k } are linear combinations of correctable 
errors with respect to the code Pq. Since each coupling operator A k acts only on 0(1) 
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qubits, the distance of the code P need to be larger than some constant value depending 
on the locality of the coupling Hamiltonian H int . Thus formally the bound Eq. fl68l) 
applies even to microscopic systems that consist only of a few qubits. Note however that 
the degeneracy of the ground-space for microscopic systems is not stable under small 
perturbations of the memory Hamiltonian H. Making the ground-space degeneracy 
insensitive to perturbations requires codes with a macroscopic distance which can be 
achieved only for macroscopic systems. 
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